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Abstract 

In previous papers, we have proposed a model that the high-frequency quasi- 
periodic oscillations (HF QPOs) observed in black-hole and neutron-star X-ray bi- 
naries are inertial-acoustic oscillations that are resonantly excited on a one-armed 
deformed disk by nonlinear couplings between the oscillations and the disk deforma- 
tion. In this paper we show that in addition to the inertial-acoustic waves, one-armed 
corrugation waves are also excited in the deformed disks. They are low-frequency os- 
cillations. We examine frequencies and their correlations among the inertial-acoustic 
oscillations and the corrugation waves that are excited, in order to know whether they 
can describe observed frequency correlations among kHz QPOs and low-frequency 
QPOs (LF QPOs) in neutron-star X-ray binaries. The results seem to well describe 
the observed correlations in Cir X-l, if we adopt M = 1.5 ~ 2.0M Q and a* ~ 0.8, 
where M is the mass of the central star and a* is the dimensionless spin parameter 
of the metric. Finally, assumptions involved in this disk-oscillation model are briefly 
summarized and discussed. 
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1. Introduction 

In previous papers (Kato 2004, 2008), we showed that in a deformed relativistic Keplerian 
disk, inertial-acoustic oscillations (and g-mode oscillations) are resonantly excited. The defor- 
mation required for the excitation is a warp or a one-armed plane-symmetric spiral pattern. 

The inertial-acoustic oscillations resonantly excited in the above deformed disks well 
describe the 3:2 twin high-frequency quasi-periodic oscillations (HF QPOs) observed in black 
hole X-ray binaries, if the disk region where the oscillations are excited is surrounded by a 
hot corona (torus) (Kato and Fukue 2006). In kHz QPOs in neutron-star X-ray binaries, 
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the frequencies of the twin QPOs change with time with correlation. These correlated time 
change of twin kHz QPOs can be accounted for, if the deformation has a time-dependent 
precession (e.g., Kato 2007)0. Any successful models of QPOs of neutron-star X -ray binaries, 
however, should account for the observed frequency correlations not only between twin kHz 
QPOs but also those among kHz QPOs and low-frequency QPOs (LF QPOs) (Psaltis et al. 
1999; Boutloukos et al. 2006). 

Concerning disk oscillation modes that are excited in one-armed deformed disks, we find 
by a careful inspection of mathematical stability analyses by Kato (2008) that the one-armed 
c-mode oscillations are also resonantly excited in addition of inertial-acoustic oscillations (and 
g-mode oscillations) (see subsection 4.2). This is of interest in relation to observations, since the 
one-armed c-mode oscillations have low frequencies and thus they are one of possible candidates 
of low-frequency QPOs (LF QPOs). Disk oscillations that are excited in one-armed deformed 
disk are thus inertial-acoustic oscillations (and g-mode oscillation) and c-mode oscillations. 

In this paper, we examine frequency-frequency correlations among the disk oscillations 
that are resonantly excited, in order to compare them with observed correlations of QPOs. 
The oscillation modes considered here are three oscillation modes, i.e., two basic inertial- 
acoustic oscillations among excited inertial-acoustic ones and the one-armed c-mode oscillation. 
Comparisons with observations are made to Cir X-l, since frequencies of twin kHz QPOs and 
LF QPOs to this source are carefully examined recently (Boutloukos et al. 2006). 

2. Overview of Resonantly-Excited Disk-Oscillation Model 

We present here the essential part of the model (Kato 2004, 2008), although the model 
has various variations. The disk we consider is a geometrically thin, relativistic disk. 

2.1. Disk Deformation 

The disks are assumed to be deformed from an axially-symmetric steady state by some 
external or internal cause. The deformation is taken to be eccentric in the disk plane, i.e., 
ma = 1, where rrid is the azimuthal wavenumber of the deformation. Concerning the symmetry 
of the deformation with respect to the equatorial plane, two cases are considered. One is the 
case where the deformation has no vertical asymmetry, i.e., n& = 0, where n<j is zero or a positive 
integer specifying the number of node(s) in the vertical direction. The other is the case where 
the deformation is asymmetric with respect to the equator, i.e., n& = \. The latter deformation 
(i.e., rrid = 1 and = 1) is a kind of tilt or warp. 

The deformation is assumed to have a time-dependent precession, i.e., it rotates in 
the azimuthal direction with time-dependent angular velocity. The angular velocity of the 
precession is denoted by u p . If uj p is positive the precession is prograde, while it is retrograde 

1 In the paper by Kato (2007), mode identification of disk oscillations to observed QPOs was not relevant. 
This will be corrected in this paper. 
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when Up < 0. In the case of neutron star, unlike the case of black hole, it has a surface, and 
the disk surrounding the star has strong radiative and magnetic couplings with the star. We 
assume that this is one of possible causes of time- dependent precession of disks. For example, 
Pringle (1992, 1996) and Maloney et al. (1996) have shown that radiative force from the central 
star can lead to corrugation of the disk. Even initially planar disks are unstable to warping 
(Pringle 1996). Movies made by Pringle and his collaborators show that the disk precession is 
time-dependent and it can drastically change even the direction of the precession. 

2.2. Disk Oscillations 

We consider oscillations on the disks described above. Disk oscillations are generally 
described by (u, m, n), where uj and m are the angular velocity and azimuthal wavenumber 
(m = 0,1,2,...) of the oscillations, respectively, and n is an integer {n = 0,1,2,...) describing 
the number of node(s) of the oscillations in the vertical direction (see Kato et al 2008 for 
detailed classification of disk oscillations). For a set of (a;, m, n), there are two different kinds 
of modes of oscillations, except for the case of n — 0. In the case of n — 0, we have inertial- 
acoustic oscillations (p-mode) alone, which consist of nearly horizontal motions (motions in the 
disk plane). In the case of n > 1, we have two different modes of oscillations. One is gravity 
oscillations (g-mode). The other is corrugation waves (c-mode) (n = 1) or vertical-acoustic 
oscillations (vertical p-mode) (n > 2). 

The oscillation modes that are treated in this paper are the inertial-acoustic mode (n = 0) 
[and the g-mode § (n > 1)] and the c-mode oscillations [n = 1) with m = 1. 

2. 3. Nonlinear Resonant Processes 

A nonlinear coupling between the disk deformation characterized by (co> p , m^, n^) and an 
oscillation with (u, m, n) brings about oscillations described by (lj±lj p , mim^, nin^), where 
various combinations of ± are possible. These oscillations are hereafter called intermediate 
oscillations. These intermediate oscillations have resonant interaction with the disk at particular 
radii where the dispersion relation of the intermediate oscillations is satisfied. There are two 
kinds of resonance for the same set of (u±lj p , m±md, n±nd). In one of them, resonance occurs 
through horizontal motions, while it occurs through vertical motions in the other one. We call 
the former resonance a horizontal resonance, while the latter resonance a vertical resonance 
(Kato 2004, 2008). 

After making the resonant interaction with the disk, the intermediate oscillations non- 
linearly couple with the deformed part of the disks to feedback to the original oscillation of (u, 
m, n) (see figure 1 of Kato 2004). This nonlinear feedback processes amplify or dampen the 
original oscillations, depending on the types of oscillations and of resonance. Careful stability 



2 In this paper we implicitly consider the g-mode oscillations together with the inertial-acoustic oscillations, 
since they are treated in the same mathematical framework as the inertial-acoustic waves, except when we 
consider their propagation region. 



analyses have been done by Kato (2008) (see also Kato 2004) in the case where the deformation 
has no precession, i.e., u p = 0. 



3. Conditions of Resonance and Excitation 

Here, we outline how the results derived by Kato (2008) to the case of uj p = are 
generalized to the case of uj p ^ 0. 

3.1. The radius of resonance 

As mentioned in the previous section, there are two kinds of resonance, i.e., horizontal 
and vertical resonances. Hereafter, we restrict our attention only to the horizontal resonance, 
since in the case of no precession, the resonance that can excite disk oscillations is the horizontal 
resonance alone (Kato 2004, 2008). In the case of uo p = 0, the condition of the horizontal 
resonance is [u — (m ± l)fi] 2 — k 2 = (Kato 2004, 2008), where Q is the angular velocity of 
disk rotation and taken to be the relativistic Keplerian one, when its numerical figure is 
required, since we consider geometrically thin disks. The symbol k is the epicyclic frequency. 
If the deformation has a precession, the condition for horizontal resonance is generalized to 



where arbitrary combination of ± is allowed. Equation (TjQ) gives resonant radii as functions of 



To find whether the resonance excites or dampens the oscillations, we must examine (i) 
the sign of work done on the oscillations at the resonant radii and (ii) the sign of the wave 
energy of the oscillations (Kato 2008, see also Kato 2004). For example, if a positive work is 
done at a resonant radius on an oscillation of a positive energy, the oscillation is amplified. 

3.2. Work done on oscillations 

Analyses made by Kato (2004, 2008) show that the sign of the work done by resonance 
on oscillations at a resonant radius is governed by the sign of u — (m ± 1)Q at the resonant 
radius, when the deformation has no precession. Let us denote resonant radii given by equation 
CQ) by r rcsonance and the work done on oscillations by the resonance at r resonancc by W. Then, in 
the case where the deformation has a precession with angular frequency u p , we have 



where sign(A) represents the sign of A, and the subscript resonance denotes the value at 

'"resonance • 

The wave energy, E, of oscillations depends on whether the main part of oscillations 
exist inside or outside the corotation radius defined by w- mfi = 0. If it exists inside the 
corotation radius, the wave energy is negative. In the case where E > 0, the oscillations are 
excited by resonance, if the work done on the oscillations by resonance is positive (W > 0). In 




(1) 



uj for a given set of u p and m (see figures 1 and 2). 




(2) 
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the case of E < 0, on the other hand, W < excites the oscillations. 

3.3. Propagation Region of Oscillations 

For oscillations to be excited efficiently, the resonant radius must be in the propaga- 
tion region of the oscillations. In the case of inertial-acoustic oscillations with uj and m, the 
propagation region is described by 

(u-mtt) 2 -K 2 >0, (3) 

while it is 

(u - mtt) 2 - k 2 < (4) 

in the case of g-mode oscillations. 

In addition to the inertial-acoustic oscillations and g-mode oscillations, we consider in 
this paper the c-mode oscillations (n—1) with m — 1. Their propagation region is described 
by (e.g., Kato 2001, Kato et al. 1998, 2008) 

(u-tt) 2 -tt 2 ± <0, (5) 

when local approximations are adopted, where fl± is the vertical epicyclic frequency (Kato 
1990) and slightly smaller than Q(= Qk)- 

4. Description of Resonantly-Excited Oscillations on Propagation Diagram 

Descriptions in the previous section give us necessary materials for specifying the os- 
cillation modes that are really excited by resonance. We now examine this issue by using 
the so-called propagation diagram of oscillations. If an arbitrarily large precession is allowed, 
various oscillations can be excited. In order to simplify situations by excluding less realistic 
cases, we restrict our attention here only to the oscillations that can be excited even in the 
limit of no precession, and study how the frequencies of such resonantly-excited oscillations are 
affected by the presence of precession of deformation. That is, we restrict our attention only 
to horizontal resonance. C^es of inertial-acoustic osculations (and g - mo de oscillation^ and 
c-mode oscillations are considered separately below, since the resonant radii are different in 
these cases. 

4-1. Inertial-acoustic oscillations 

Examination of excitation conditions of inertial-acoustic oscillations on the propagation 
diagram has been made in a previous paper (Kato 2008) in the case of no precession (see figures 

3 The resonant excitation of the g-mode oscillations can be treated together with that of the inertial-acoustic 
oscillations, as mentioned before, except for the difference of propagation region. Excitation of g-mode 
oscillations is, however, less interesting compared with that of the inertial-acoustic oscillations, since in their 
propagation region the co-rotation resonance appears except for special cases, and it dampens the g-mode 
oscillations (Kato 2003; Li ct al. 2003). Hence, hereafter, we do not discuss the g-mode oscillations. 
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Fig. 1. Propagation diagram (radius- frequency relation) showing the radii and the frequencies of iner- 
tial-acoustic oscillations (with m = 1) that resonantly interact with a deformed relativistic Keplerian disk. 
The disk deformation required is a warp or a one-armed plane-symmetric deformation. The mass of the 
central source is taken to be 2.4M© with no spin. Two case are shown where the deformation has no 
precession (uj p = 0) and has a precession of lo p = —100 Hz. The inside area bounded by two curves labeled 
by il + K and Q — k are the evanescent region of the oscillations. The thick curves (there are two curves 
for each w p , and in the case of Ct> p < they further appear in the right-lower corner of the figure) show the 
radius-frequency relation of the oscillations that arc excited. The thick curves continue in the evanescent 
regions of oscillations. In the evanescent region, however, excitation of the oscillations will not occur 
efficiently, and hence the curves in the region are shown by thick but dotted curves. At the boundaries 
between the propagation and evanescent regions the excitation will occur most efficiently, since the oscilla- 
tions stay there for a long time because of vanishing of the group velocity. The thin dotted curves show the 
radius-frequency relation of resonance where oscillations are damped. In the case of gravity oscillations, 
the propagation region and the evanescent region are opposite to those of the inertial-acoustic oscillations. 
That is, the evanescent region of the inertial-acoustic oscillations becomes now the propagation region. In 
other words, in the case of gravity oscillations, the thick dotted curves in figure 1 should be changed to 
thick curves, and the thick curves to thick dotted ones. 

1 and 2 of Kato 2008). Here we extend it to cases where the deformation has precession. 

Figure 1 is the frequency-radius relation in the case where oscillations are one-armed 
(i.e., m = 1). The central source is taken to have no spin (a* = 0) and the mass M is 2.4M , 
where a* is the dimensionless spin parameter specifying the rotation of the metric and a* = 1 
corresponds to the extreme Kerr. The inner region bounded by two curves labeled by uo = Q + n 
and uj = Q — k is the evanescent region of the one-armed inertial-acoustic oscillations [cf., 
inequality (J3J)], and the outside of the region is the propagation region of the oscillations. Thick 
(both solid and dotted) curves give the frequency-radius relation for resonant excitation of 
one-armed (m = 1) inertial-acoustic oscillations [see equation (CQ) and examine the signs of W 
and E\. The dotted parts of the thick curves are, however, inside the evanescent region of the 
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oscillations [see equation and thus excitation on the part will be practically inefficient. 
That is, the frequency-radius relations for excitation are the thick solid curves. Two cases of 
uj p = and u p = — 100 Hz are shown. The thin dotted curves are the frequency-radius relation 
of resonance that leads the oscillations to damping. 

As shown in figure 1, one-armed oscillations that are excited are not unique. In a finite 
range of radius, oscillations are excited with different frequencies [see the sentence just below 
equation flTJ)]. Furthermore, two oscillations can be excited at the same radius. One has higher 
frequencies than Q at the resonant radius, while the other is smaller than Q. 

Here, we assume that the oscillations whose group velocity is slower will be excited till a 
larger amplitude than those with a faster group velocity, since they stay a place for a longer time 
and the resonance condition is maintained for a longer time. If we adopt local approximations, 
the group velocity of oscillations vanishes at the boundary between the propagation region and 
the evanescent region, which is specified by uj = Q + k or u = Q — k in the case of one-armed 
oscillations. If this consideration is adopted, the most predominantly excited oscillations are 
those with uj = Q + k in the case of the higher frequency oscillations and those with uj = Q — k in 
the case of the lower frequency oscillations. Then, combining this condition and the resonant 
condition given by equation flTJ), we have, as the condition determining the resonant radius, 

K = -(n-u p ). (6) 

If the deformation has no precession, uj p = 0, this condition gives k = fi/2 and the resonance 
occurs at 4r g when the central source has no spin, where r g is the Schwarzschild radius given 
by r g = 2GM/c 2 . Among a series of resonant radii given by equation (JTJ, the resonant radius 
given by equation (J6j) is particularly denoted hereafter by r res , which is a function of a*, M, 
and Up. 

If we use the notations used so far in the series of our papers, the lower-frequency 
oscillation at r rcs is cull an d the higher-frequency one is co>h, i.e., 

u LL = (Q - k) tcs and uu = (Q + K) ieB . (7) 

Next, we proceed to the two-armed oscillations (m = 2). Figure 2 is the same as figure 
1 except for m = 2. Two cases of uj p = and uj = —100 Hz are again shown for M = 2AM & 
and a* = 0. The thick solid curves and the thick dotted curves are the frequency-radius re- 
lations for resonant excitation of oscillations. The thick dotted parts are, however, inside the 
evanescent region of the oscillations and thus the oscillations in that part will be practically 
not excited. Thin dotted curves represent the resonance that dampens the oscillations. As 
in the case of oscillations of m = 1, the most prominently excited oscillations will be those in 
which resonance occurs at the boundary between the propagation region and the evanescent 
region. This consideration again leads to equation ([6]). Following again the notation used in 
the previous papers, we denote the frequency of the lower-frequency oscillations at r res and that 
of the higher- frequency oscillations by u;l and cuhh, respectively, i.e., 
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Fig. 2. Same as figure 1 except for m = 2. In the present case of m = 2, the evanescent region of the 
oscillations is the inside region bounded by two curves of 2Q + n and 2£1 — k (these curves are not labeled 
in the figure in order to avoid complexity). 



As figures 1 and 2 show, the curves of the radius-frequency relations for excitation of 
oscillations (i.e., thick solid or thick dotted curves) cross again, in an outer radius, the boundary 
curve between the propagation region and the evanescent region when uj p = —100 Hz. In order 
to show this, the resonant radius, r res , obtained by solving equation (jHJ) is shown in figure 3 as 
a function of u> p for a* = 0.1 and a* = 0.3 in the case of M — 2.OM . In figure 3 the r Tes -u p 
relation for the corrugation waves of m — 1 is also shown (see the next subsection). [In the 
case of corrugation waves, the resonant radius is different from r res . That is, it is not given by 
equation (jfJJ), but by equation (Q (see the next subsection).] 

4-2. C-mode Oscillations 

Here, we consider excitation of one-armed corrugation waves (c-mode oscillations) re- 
sulting from horizontal resonance. The c-mode oscillations that we consider here are those 
of m = 1. They are a kind of tilt (warp) and have very low frequencies. Their propagation 
region is specified by equation ([5]). The boundary of the propagation region, to = Q — Q±, in 
the frequency-radius diagram is shown in figure 4 by a thin solid curve for two cases of a* = 0.3 
and a* = 0.8 with M = 2.OM . In figure 4, the left-hand side of the curve of uj = Q — Q± is 
the propagation region of the one-armed corrugation waves. The right-hand side of the curve 
is the evanescent region . 

The horizontal resonance occurs at the radii where [(u ± u p ) — (m ± l)fi] 2 — k 2 = is 
satisfied as mentioned before [see equation (CQ)]. Examination of the sign of work done on the 




and 




(8) 
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Fig. 3. The relation between resonant radius (r rcs or r corr ) and uj p . The relation is shown for inertial-a- 
coustic oscillations and the c-mode oscillations with m = 1 . Two cases of spin parameter of the central 
source, i.e., a* = 0.1 and a» = 0.3, are shown with M = 2.OM0. It is noted that there are two resonant 
radii for a given uj p in the case of uj p < 0. For the case of a* = 0, see figure 1 by Kato and Fukue (2007). 
In the case of a* = 0, the curve corresponding to the corrugation waves with m = 1 disappears in the 
approximation used in this paper. 




Fig. 4. Propagation diagram showing the radius - frequency relation for resonantly excited corrugation 
waves with m=l. Two cases of a* = 0.3 and a* = 0.8 are shown for two cases of precession, lu p = and -100 
Hz, with M = 2.0M©. The curves labelled by il — fl± are the boundary between the propagation region 
and the evanescent region. The left-hand side of the curves are the propagation region. As in figures 1 and 
2, the thick solid curves represent the radius - frequency relation where oscillations arc excited. Although 
not shown in this figure, the thick solid curves appear again in the outer region of disks when uj p = —100 
Hz (see figure 5). 
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Fig. 5. Enlargement of the lower-right corner of figure 4. Notice that in the case of lu p = —100 Hz, the 
dotted curve crosses the curve of lu = il — f2j_. Outside the crossing point, the dotted curve changes into 
a thick solid curve (i.e., resonance excites oscillations). 

one-armed (m = 1) corrugation waves at the resonant radii given by equation ([1]) and the sign 
of the wave energy of the waves shows that among resonances specified by equation ([I]) the 
resonance given by uj = k + uj p excites the oscillations. In figure 4, the curve of u = k + uj p 
is denoted by a thick solid curve when it is in the propagation region (uj < Q — Q_L) of the 
oscillations, while by a dotted curve when it is in the evanescent region. Although two cases of 
uj p = and uj p = —100 Hz are shown in figure 4, the difference of two curves of uj = k + uj p in 
two cases of uj p = and uj p = —100 Hz is little, when the curves are in the propagation region 
of the oscillations. In the right-lower corner of figure 4, the curve of uj = k + uj v crosses again 
the curve of uj = Q — Q± when uj p < 0, and the curve of u = k + u p enters into the propagation 
region of the corrugation waves for large value of r. In order to show this, an enlargement of 
the right-lower corner of figure 4 is shown in figure 5. 

As in the case of inertial-acoustic oscillations, the group velocity of the c-mode oscilla- 
tions vanishes at the boundary of the propagation and evanescent regions. Hence, we consider 
that the c-mode oscillations predominantly excited by the resonance are those that satisfy 
simultaneously uj = Q — Q± and uj = n + uj p . Combination of these two equations gives 



This is the resonant condition for the most predominantly excited c-mode oscillations. The res- 
onant radius determined by this condition is denoted by r corr . The frequency of the corrugation 
waves at this resonant radius is denoted by uj covt , i.e., 



Equation (Q is satisfied at two radii when uj p < as is shown in figures 4 and 5. The radius- 
precession relation specified by equation ([9]) is shown in figure 3 in the case of a* = 0.1 and 0.3 




(9) 




(10) 



with M = 2.0M 



10 



5. Frequency Correlations and Comparison with Observations 

Arguments in subsection 4.1 show that the inertial-acoustic oscillations that are dom- 
inantly excited on the deformed disks are those whose frequencies are specified by cull, ^l, 
and wh [see equations ([7j) and flSJ)]. (Oscillations with u; H h are not considered here, since their 
frequency is rather high and will be outside of observations.) Among these three oscillations, 
the oscillation of c^h will be less important in observational view points. This is because the 
propagation region of the oscillation is outside the resonant radius r res and extended infinitely. 
This suggests that the oscillation will not remain for a long time around the resonant radius, 
since it propagates away outside. On the other hand, the propagation regions of waves with 
co>ll and c^l are inside the resonant radius. The waves are, thus, partially trapped between the 
resonant radius, r res , and the inner edge of the disk, and remain in that region for a long time, 
different from the wave with un. 

Next, we should emphasize here that one-armed oscillations with cull are not always 
observed only with their own frequency c^ll- They may be observed with the two- fold frequency 
when the disk region where the oscillations are excited is surrounded by a hot corona and when 
we observe high energy photons that are reprocessed in the corona (see figures 2 - 4 by Kato 
and Fukue 2006). Based on these two considerations, we think that the frequencies related to 
the twin HF QPOs observed in black-hole and neutron-star binaries are 2u;ll and u^. 

The corrugation waves that are excited are also one-armed, and thus by the same reason 
mentioned above they will be observed with the twofold frequency, 2u corr . It is further noted 
that they are trapped inside the resonant radius, r corr . 

In summary, we take the picture that the main frequencies to be observed are 

2cj corr , 2u) LL , (11) 

in the order of increase of frequency. As discussed below, 2u co „ will correspond to the low- 
frequency QPO (LF QPO), and 2cjll and c^l correspond to the lower and upper kHz QPOs, 
respectively. 

If the deformation has no precession, the frequency ratio of 2u;ll to cul(= wh) is just 2 
to 3, independent of M and a*. We think that the twin QPOs in black-hole X-ray binaries 
represent this case (Kato and Fukue 2006). 

In the case of neutron-star X-ray binaries, we suppose that the disk deformation has 
a time- dependent precession as mentioned before (see subsection 2.1). The frequencies 2u corr , 
2o;ll, and are thus functions of u p as well as M and a*. Then, rearranging the functional 
relations, we obtain 2u COTT , u^, and \u p \ as functions of 2o>ll with parameters a* and M. As a 
typical example, these relations are shown in figure 6 for M = 2.0M© and a* = 0.3. It is noted 
that this figure is drawn for uj p < 0. When u p < 0, we have two resonant radii for a given uj p 
(see figure 3), i.e., an inner resonance and an outer resonance. The curves drawn in figure 6 
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Fig. 6. Diagram showing the wl - 2wll and 2u> corl - - 2wll relations in the case of M = 2.0Mq and 
a* = 0.3. The curve of the |w p | - 2ull relation is also shown by thin curve. Observational data showing 
the correlation between the upper and lower kHz QPOs and the correlation between the low frequency 
QPO and the lower kHz QPO (Boutloukos ct al. 2006) are superposed. The filled triangles are the data 
of Cir X-l. The open circles are data of other sources of Psaltis et al. (1999), which have been plotted on 
Boutloukos ct al.'s figure. 
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Fig. 7. The same as figure 6, except that two cases of (M = 1.5M , a* = 0.8) and (M = 2.0M Q , a* = 0.8) 
are drawn on the same diagram. The blue curves are for (M = 1.5M Q , a» = 0.8) and the red ones are for 
(Al = 2. OA/©, a» = 0.8). The green curve represents the precession frequency. The differences between the 
two curves of 2w corr are a few, except in the high frequency region. 
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Fig. 8. Diagram showing the Aw ~wl relation, where Aw is the difference between wl and 2wll, i.e., 
Alo = lj-l — 2wll = ^rcs- Parameters (M, a») adopted are, from the uppermost to lowermost curves, (2.0, 
0.8), (2.0, 0.3), (2.0, 0.0), and (2.4, 0.0). The diagram showing the difference between the upper and lower 
kHz QPO frequencies versus the upper kHz QPO frequency (Boutloukos et al. 2006) are superposed, 
assuming that the upper kHz QPO corresponds to ui-^. The black circles (they are on the left-hand side 
of this figure) are the data of Cir X-l (Boutloukos ct al. 2006) and the dark circles are the data of other 
sources (Psaltis et al. 1999). The thin solid curve is the relativistic precession model by Stella and Vietri 
(1999) (sec Boutloukos ct al. 2007). 

(and in subsequent figures) are for the frequencies of the outer resonance^. On figure 6 the 
observational data of Cir X-l (figure 9 of Boutloukos et al. 2006) are superposed, assuming 
that the frequency of the lower kHz QPO corresponds to 2co>ll- The frequency correlations 
shown in figure 6 seem to qualitatively account for the correlations among observed QPOs of 
Cir X-l, but there are systematic deviations, especially in the correlation between the 2u; corr - 
2o;ll relation and the observed LF QPO - lower kHz QPO relation. We have two adjustable 
parameters, i.e., M and a*. The 2ui COTT - 2c<jll relation depends strongly on the value of a*, since 
uj corr is related to a small difference between two large quantities of Q and Q±, and Q± becomes 
close to Q with decrease of a*. On the other hand, the u;l _ 2o;ll relation depends only weakly 
both on a* and M . The observed data of Cir X-l seem to be well described, if M = 1.5 ~ 2.OM 
and a* = 0.8 are adopted, as shown in figure 7. It is noted here that if we adopt a* = 0.6 with 
M = 1.5 ~ 2.OM , the 2uj co „ - 2o;ll relation runs slightly below the observed LF QPO - lower 
kHz QPO relation. 

In order to see more closely in what parameter cases the uo^ _ 2c<jll relation fits the 
observations, the difference, Au;, between o>l and 2u;ll, i-e., Acu = lu l — 2uj ll = /t res , is plotted 
against u;l in figure 8 for M = 2.0M Q with some values of a*. For a comparison, the case of 
M = 2AM & and a* = is also shown in figure 8. The observed frequency difference between 

4 The inner resonance has higher frequencies and is inadequate to describe the QPO frequencies observed in 
Cir X-l. In some sources, kHz QPOs have been observed with higher frequencies, compared with in Cir X-l. 
Such high frequency QPOs will come from the inner resonance, as discussed in the next section. 
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Fig. 9. The same as figure 8 except that M = 1.7 M® and a* = 0.8, 0.3, 0.0 from upper to lower. 

the upper and lower kHz QPOs versus the upper kHz QPOs in Cir X-l is given in figure 11 
of Boutloukos et al. (2006). The figure of Boutloukos et al. (2006) has been overlapped on 
figure 8. As another example of comparison, the cases of M = 1.7 Mq are shown in figure 9 for 
three cases of a* = 0.8, 0.3, and 0.0. Figures 8 and 9 show that a rather high value of a*, say 
a* ~ 0.8, is necessary in the case of Cir X-l in order to account for the correlation between the 
twin kHz QPOs. Requirement of such a high value of a* is compatible with what required in 
the previous paragraph to describe the LF QPO - lower kHz QPO correlation. 

6. Discussion 

We have applied a resonantly-excited disk-oscillation model to Cir X-l in section 5. A 
comparison of observational data with the model suggests that the spin of the central source of 
Cir X-l is rather high, i.e., a* ~ 0.8. Let us roughly estimate the spin frequency u B of Cir X-l by 
using the estimated value of a*. If a neutron star is rotating uniformly inside the star with u s , 
the specific angular momentum, a, of the star is given by a = 2nu s I/cM, where / is the moment 
of inertia of the star. Since a* is related to a by a* = (c 2 /GM)a, we have u B = a*(GM 2 /2%cl), 
which gives 



where M = M/M Q and J 45 = J/10 45 . If we take M = 2, M // 4 5 = 1 (e.g., Stella and Vietri 
1999), and a* = 0.8, we have u s ~ 2.2 x 10 3 Hz. Such a high spin, as well as a high eccentric 
orbit of Cir X-l, may suggest that the system is young. 

Comparison of observations and the model suggests that the disk of Cir X-l has a 
retrograde precession and it varies in the range of 100 Hz ~ 200 Hz (see figures 6 and 7). For 
kHz QPOs to be excited by retrograde precession, the absolute value of the precession must be 
less than a few hundred Hz (see figure 3). Otherwise, inertial-acoustic oscillations cannot have 




(Hz), 



(12) 
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resonant interaction with disk deformation, i.e., we have no kHz QPOs. That is, the range 
of variation of QPO frequencies is related to the magnitude of precession and the range of its 
variation. 

Here, we more closely discuss the relations among precession, resonant radius, and fre- 
quency of resonant oscillations. If precession is retrograde and its absolute value is smaller 
than a critical value, resonance occurs at two different radii for each of inertial-acoustic and 
corrugation waves (see figure 3). For each case of inertial-acoustic and corrugation waves, we 
call the resonance at the outer radius an outer resonance, and the resonance at the inner radius 
an inner resonance (see section 5). In the case of Cir X-l, as shown in section 5, the resonance 
that can describe observations is the outer one (see also that the gradients of three curves of 
\uj p \, cul, and 2lj cott are all positive in figures 6 and 7). In Cir X-l, the frequencies of oscillations 
resonantly excited at the inner resonance are higher than those observed, as long as conven- 
tional mass of neutron stars are adopted. The fact that oscillations at the inner resonance are 
not observed in Cir X-l suggests that the inner disk region of Cir X-l is highly disturbed. One 
of possible causes of such disturbances will be strong couplings of the inner disk region with 
the surface of the central source through strong magnetic fields and high spin of the central 
source. 

The cases where the observed twin kHz QPOs are the oscillations excited at the outer 
resonance will be rare, compared with the cases where they are those excited at the inner 
resonance, since the twin QPOs in many sources have higher frequencies than those in Cir 
X-l. That is, in the sources where the high-frequency twin QPOs have really frequencies close 
to or higher than one kHz, they will be resonant oscillations at the inner resonance and the 
precession will be small or prograde. It is noted that in the case where the precession is slow 
or prograde, the outer resonance will be observationally less interesting, since the frequencies 
excited there are too low or the outer resonance itself is absent (see figure 3). 

Let us briefly demonstrate an example in which QPOs occur at the inner resonance. Sco 
X-l has twin kHz QPOs in the 800 - 1100 kHz range, a ~45 Hz horizontal-branch oscillation 
(HBO), and a ~ 6 Hz normal-branch oscillation (NBO). Here, we examine the possibility that 
the observed twin kHz QPOs correspond to 2cjll an d and the 45 Hz oscillation to 2u corr , 
as in the case of Cir X-l. The observed frequencies and their changes of the lower and upper 
kHz QPOs can be roughly described by assuming that the precession of the disk deformation 
is prograde and M ~ 2.4M , when a* is taken to be zero (see figure 3 of Kato 2007). If a* is 
zero, however, there is no corrugation wave (uj corr = 0) and we cannot account for a ~ 45 Hz 
oscillation by our present model. As mentioned in section 5, the frequency correlation between 
2uj ll and is little affected by changing the value of a*, while the frequency of 2u corr is 
strongly affected by a*. Hence, by adjusting the value of a* we can describe a 45 Hz QPO by 
2uj corr . Our preliminary results show that both frequencies of twin kHz QPOs and of HBO can 
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be accounted for by our model with M ~ 2.6M(jj and a* ~ 0.15, the precession being nearly 
zero or weakly prograde. A ~ 6Hz normal-branch oscillation (NBO) might be a manifestation 
of the slowly precessing disk deformation required in our model. 

It is noted here that we have adopted 2u; corr (not u corr ) and 2o>ll (not cull) as frequencies 
of LF QPO and the lower kHz QPO, respectively. This is, however, not always the case in all 
sources, unless the disks are surrounded by coronae. Hence, some caution will be necessary 
when we make comparison of the model with observations. For example, the frequency of LF 
QPO may be u; C orr in some objects. 

Hereafter, we mention the basic assumptions involved in the model, and discuss validity 
of these assumptions. First, we have assumed that (i) the disk is deformed in an eccentric form 
(ma = 1) in the state when the QPOs are observed. There seems no clear observational evidence 
for or against this assumption. There are, however, numerical simulations that support a rela- 
tion between disk deformation and appearance of QPOs. Magnetohydrodynamical simulations 
of accretion flows by Machida and Matsumoto (2008) show a formation of a torus in the central 
region of disks and its quasi-periodic deformation into a one-armed form at a certain phase of 
the flows. They show that this deformation gives rise to a low-frequency quasi-periodic varia- 
tion of mass accretion rate, and further that in this deformed state, high-frequency QPOs of 
the order of one hundred Hz appear (the mass is taken to be 10 M & ). It is not clear, however, 
whether such disk deformations are related to the disk deformation required in the present 
disk-oscillation model. 

The second requirement involved in the present model is that (ii) the deformation must 
have a time-dependent precession. This time-dependent precession is necessary to bring about 
a change of resonant radius and thus to account for the frequency change of QPOs. It is 
uncertain whether such precession really exists in neutron-star binary sources. By couplings 
through magnetic and radiation fields, the accretion disks surrounding a neutron star will 
have strong interaction with the surface of the star. This may be one of possible causes of 
time-dependent precession of deformed disks (see subsection 2.1). The Papaloizou - Pringle 
instability in a torus, however, will not be at least the cause of the disk precession in Cir X-l, 
since the precession required to describe the observations is retrograde. The hectohertz QPOs 
might be a manifestation of the precession, but if so, it is not clear why they are not always 
observed in the phases where QPOs appear. 

Third, it is assumed that (iii) one-armed (m = 1) disk oscillations excited on disks are 
observed with the two-fold frequency. That is, we assume that the resonant region of the 
geometrically thin disks where QPOs are excited is surrounded by a hot corona (a torus) and 
the observed QPO photons are those reprocessed in the corona to higher energy. If such a 
picture is taken, we can expect the two-fold frequency in the case of one-armed oscillations (see 

5 This value of M is rather high compared with the conventional one adopted as the mass of neutron stars, 
but will be still in the range allowed theoretically and observationally. 
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figures 2 - 4 by Kato and Fukue 2006). At the phase when the QPOs are observed in black hole 
sources (i.e., the very high state), the spectra really show coexistence of a disk component and 
a corona (a torus). It is uncertain, however, that such corona (torus) exists even in the case of 
neutron-star X-ray sources, since in the case of neutron stars, a hot corona will be weakened 
by presence of abundant soft photons from the surface of the star. 

Fourth, an important assumption involved in our model is that (iv) dominantly excited 
oscillations are those with zero group velocity. This is an approximate procedure. In the 
future, it will be necessary to solve global disk oscillations with relevant boundary conditions 
in order to clarify whether this approximate procedure is allowed as the first approximation. 
Another important issue to be reminded here is whether the outward propagation oscillations, 
characterized by c^h and u>hh, can be really regarded as less important, compared with the 
oscillations characterized by c^ll and u-^. 

Unlike the epicyclic resonant model by Abramowicz and Kluzniak (2001) and Kluzniak 
and Abramowicz (2001), in the present disk-oscillation model, the appearance of the twin 
QPOs is not a result of mutual interaction of twin oscillations. In the present model, twin 
oscillations are independently excited by interaction with the disk deformation. In this reason, 
there is no direct amplitude relation between the twin oscillations in the present model. Torok 
(2008) found a frequency - amplitude relation of the observed twin QPOs. Such correlations 
cannot be accounted for in this model, unless nonlinear evolutions of the excited oscillations 
are considered. 

Finally, we should notice a close relation between our present model and the relativistic 
precession model by Stella and Vietri (1999). In the latter model, the upper and lower kHz 
QPOs are interpreted, respectively, as a frequency of disk rotation, Q, and periastron precession 
frequency, Q — k, at a certain radius, which is a free parameter. In the present disk-oscillation 
model, the frequency of the upper kHz QPO is 2Q — k (which is close to Q), and that of 
the lower one is 2(f2 — n). In the present mode, the radius specifying these frequencies is the 
radius of resonance, r res , which is determined by coupling between the oscillations and the disk 
deformation. Further, in the model by Stella and Vietri (1999), the frequency of low-frequency 
QPO is the nodal precession frequency, Q — at a radius, which is again a free parameter. 
In our model the frequency of low-frequency QPO (LF QPO) is the frequency of corrugation 
waves, which is Q — but the radius specifying the value of Q — Q± is the resonant radius 
given by r corr . 

The author thanks the referee for invaluable comments on comparison of the model 
with observations. The author also thanks M. Abramowicz, W. Kluzniak, R. Matsumoto, S. 
Mineshige, M. Bursa, J. Horak, M. Machida, G. Torok, for helpfu discussions during the YITP 
workshop YITP-W-07-14 on " Quasi- Periodic Oscillations and Time Variabilities of Accretion 
Flows". This workshop was financially supported by the Yukawa Institute of Theoretical 
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